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We examine the profile of a liquid front of a film that is 
dewetting a solid substrate. Since volume is conserved, the 
material that once covered the substrate is accumulated in a 
rim close to the three phase contact line. Theoretically, such 
a profile of a Newtonian liquid resembles an exponentially 
decaying harmonic oscillation that relaxes into the prepared 
film thickness. For the first time, we were able to observe this 
behavior experimentally. A non-Newtonian liquid - a polymer 
melt - however, behaves differently. Here, viscoelastic prop- 
erties come into play. We will demonstrate that by analyzing 
the shape of the rim profile. On a nm scale, we gain access 
to the rheology of a non-Newtonian liquid. 
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The interplay of viscous and viscoelastic properties of 
a fluid is of enormous technical importance. The rheol- 
ogy e.g. of a coating will be adjusted such that standard 
coating techniques may be applied. In those thin films, 
however, the typical rheological analysis techniques fail. 
In this Letter we will propose a method that can even be 
applied to nm thin films. We will moreover demonstrate 
the sensitivity of a liquid profile even to slight changes in 
the viscoelastic properties of the material. To the best of 
our knowledge, the analysis of a rim profile was only sub- 
ject to theoretical studies , yet experimental studies 
concentrated on the characterization of the dynamics of 
dewetting (6[jl0) , but not on the shape of the rim. 

A liquid film on a non-wettable substrate is not stable 
and will bead up: circular 'dry' regions ('holes') appear, 
the radius of which grows steadily with time . If mass 
conservation is valid, the material once at the inner part 
of the hole will be accumulated in a rim along the perime- 
ter of the hole. Close to the three phase contact line of 
the liquid, i.e. on the 'dry' side of the rim, the rim shape 
is dominated by the receeding contact angle 9 of the liq- 
uid on the surface. On the 'wet' side of the rim, where 
the actual contact angle is zero, the liquid merges some- 
how or other from the crest into the initial film thickness. 
From a theoretical point of view, the rim is a disturbance 
in the liquid and is expected to decay via an undulation, 
the wavelength of which is dominated by the surface ten- 
sion of the liquid. To the best of our knowledge, however, 
this peculiar behavior was never reported on dewetting 
experiments. In our experimental system, short chain 
polystyrene (PS) on silicon wafers, an undulatory behav- 
ior is observed and measured for the first time. It is one 
of the aims of this Letter to explain why this behavior 
is only rarely found in experiments and how it can be 



induced. Moreover, as we will show in this study, the 
profile of the 'wet' side of the rim may serve as a finger- 
print of the viscoelastic properties of the liquid. 

Using thin atactic PS films on Si wafers as a model sys- 
tem for liquid front studies has several advantages: the 
PS melt has a very low vapor pressure, hence mass con- 
servation is valid, its viscosity can be tuned by choosing 
various chain lengths and/or annealing temperatures, its 
physical properties are well-known and PS is available in 
high purity and low polydispersity. Moreover, PS films 
can be studied in the liquid state by atomic force mi- 
croscopy ( AFM) . Polished Si wafers are ideal model sub- 
strates due to their very low surface roughness (rms val- 
ues below 0.1 nm). The native oxide layer of the Si wafers 
renders the surface quite inert in air, and chemical ho- 
mogeneity of the surface can be achieved by well-known 
cleaning procedures. Moreover, to check the influence of 
contact angle on the profile of the liquid front, the wetta- 
bility of the Si wafer can be changed by the preparation 
of a self-assembled monolayer of octadecyl-trichlorosilane 
(OTS) on the surface. 

Two types of substrates were used: Type A are Si 
wafers with a thin natural oxide layer of 1.6(2) nm |l^ ] 
(Wacker Chemitronics, Burghausen, Germany; crystal 
orientation (100), p-(boron-)doped, 6-12 Vlcm) on top of 
which an OTS monolayer was prepared following stan- 
dard procedures. Type B substrates are Si wafers with a 
thick oxide layer of 191(1) nm (Silchem GmbH, Freiberg, 
Germany; (lOO)-oriented, p-(boron-)doped, conductivity 
> 1 ilcm). 

Thin polystyrene films (molecular weight between 2.05 
kg/mol and 600 kg/mol, M„/M„ = 1.02 - 1.05, purchased 
from PSS, Mainz, Germany and Polymer Labs, Church 
Stratton, UK) were spin cast from toluene solution onto 
the substrates. Before coating, type B substrates, cut in 
ca. 1 cm^ pieces, were freed from Si dust by treatment 
with a C02-jet ||l^. The samples were then 'degreased' 
by sonication in cthanol, acetone and toluene. Residual 
hydrocarbons were removed by immersing in a fresh 1:1 
II2O2 (conc.)/Il202 (30%) solution for 30 min, followed 
by a thorough rinse in hot Millipore"^*^ water. Due to 
the large contact angle of the toluene solution on type 
A substrates, PS films were first spin coated onto freshly 
cleaved mica surfaces (B & M Mica Co., Flushing, NY), 
then floated onto Milliporc"^*^ water and picked up by 
the type A substrates. 

Cleaning and coating was performed in a class 100 
clean room atmosphere. The thickness of the Si wafer's 
oxide layer and of the polymer films were measured by 
ellipsometry (Multiscope by Optrel GdBR, Berlin, Ger- 
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many) at various angles of incidence. Further character- 
ization of wafers and polymer films was done by atomic 
force microscopy (AFM) (Multimode III by Digital In- 
struments, Santa Barbara, CA) in Tapping Mode-^*^. 
Annealing took place on a temperature-controlled hot 
plate in a class 100 clean room atmosphere or in situ 
on top of the AFM sample holder. 

The inset of Fig. 1 presents an AFM scan of a typ- 
ical hole in a 6.6(2) nm ^ thick PS(2.24k) film on a 
wafer B. The diagram shows a radial cross section of this 
hole. Most outstanding feature of the cross section is the 
way the rim decays into the unperturbed film. The rim 
profile is asymmetric, with higher slopes near the three- 
phase contact line and a 'trough' on the wet side of the 
rim where it meets the undisturbed film. Here, the form 
resembles a damped harmonic oscillation. We describe 
the profile by the contact angle Q of the liquid at the 
three phase contact line, the height A of the rim at its 
crest, the depth B of the first trough and the slope tana 
of the rim profile on the wet side at the level of the initial 
film thickness Hq, as depicted in Fig. 1. The trough is 
also visible in the top view AFM scan in the form of a 
dark circle surrounding the wet side of the rim (inset of 
Fig. 1). 

From the temporal evolution of the experimental rim 
profile as shown in Fig. 2, we can determine the width a 
and height A of the rim, the depth B and the width 
b of the trough, and the slope tana as a function of 
time. Moreover, by in situ AFM scanning, we can record 
the dewetting velocity v. Within the experimental error, 
the contact angle O stays constant during hole growth, 
which is consistent with other experimental and theoret- 
ical studies |P,p^. Fig. 2b) and c) show how \B\ and a 
vary with increasing height A of the rim. From a certain 
value A*^p of the rim height onwards, which will call the 
late stage in the following, both B and a reach a plateau 
value. Below A*^p, in our terms now the early stage, \B\ 
as well as a increase - within the experimental error - 
linearly with A. 

For a physical understanding of the occurrence and the 
temporal evolution of the trough we analytically solve 
in the following the linearized thin film equation in the 
lubrication approximation [[l5|-pT[: 



(1) 



where the wavelength A is defined as 



dh ahl d^h 9h _ ^ 
dt 3ry dx"^ dx ' 



where h is the film thickness, a and rj are the surface ten- 
sion and viscosity, and v the velocity of the three phase 
contact line. We hereby neglect the disjoining pressure - 
^'(ft,), since we know from an earlier study that its impact 
on the dewetting morphology is negligible as compared to 



the interfacial forces involved here 



The 'traveling 



wave' solution /i(^) = h{x — vt) of Eq. (ll) reads 
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Here, r] denotes the viscosity of the polymer. Eq. (g) 
represents a damped harmonic oscillation. The ratio of 
/i(^) at two consecutive extrema is: 
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and is plotted as the dotted line in Fig. 3. This ra- 
tio is independent of any experimental parameter of the 
system, in particular of the dewetting velocity. It is uni- 
versal for any rim profile of a retracting Newtonian liquid 
if the film thickness is small as compared to the capillary 
length. Although, it cannot be valid for the entire pro- 
cess of dewetting. Due to accumulation of material inside 
the rim, its width a will be much larger than the width 
b of the trough. In other words, as soon as the crest of 
the rim exceeds a certain height A^j^^^ , the description of 
the entire rim profile as a damped harmonic oscillation 
with a single wavelength A=26 in Eq. (|^) will be no more 
adequate. 

In Fig. 3, the experimental data of |-B/^| are plotted 
over A. They reveal that - within the experimental error 
- B/A = const, for A < A*^^. Yet the constant is roughly 
a factor of two smaller than 0.163. For A > At 
decreases with increasing A. 

In order to extend the theoretical description to the 
late stage of dewetting, we have to introduce a boundary 
condition. From experiment, we have learned that B as 
well as a are constant in the late stage, c.f. Fig. 2b 
and c. For A > A* we therefore introduce into the 
calculation the experimental plateau value a = 1.0(1)° 
as a boundary condition. The trough depth B hence is 
simply obtained by identifying in Eq. iQ) the slope of 
/i(^) at its first zero as tana. So, B is a function of the 
experimental parameters cr, 77, ho and dewetting velocity 

V. 

Figure 3 shows the result for |-B/v4| obtained with the 
boundary condition a =1.0(1)° for A > A'^/^^^ (dashed 
line). Here, the trough depth B is calculated with a = 
30.8 mN/m, r? = 300 Pa-s (Ref. ||l|]), /iq = 11 nm, and 
v{t) taken from the online AFM measurement shown in 
Fig. 2 1^^. The point of intersection of dashed and 
dotted line at A=Al^^^ corresponds to the value of A*^^ 
at which the experimental values of are starting to 

decrease. This point hence is name A* in the following. 
For A> A* , the course of the theoretical curve for |i3/A| 
resembles that of the experimental data, but the absolute 
values differ by a factor of about 1.5. 

To sum up, one may say that the theoretical curves 
presented above well describe the qualitative develop- 
ment of the height of the rim with respect to the depth 
of the trough. A quantitative agreement, however, is 
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not reached yet. In fact, we found that the trough may 
even be completely suppressed by increasing the molec- 
ular weight of the polymer. In Fig. 4 we show profiles of 
holes in PS films on type A substrates of about the same 
film thickness (50 nm) and hole radius, yet of various 
molecular weight: 2 kg/mol, 18 kg/mol, 101 kg/mol and 
600 kg/mol. The most outstanding feature in these pro- 
files is that the trough vanishes and the width of the rim 
grows with increasing molecular weight of the polymer. 

It is interesting to note that the PS(2k) and the 
PS(18k) profiles - both polymers are below the entan- 
glement length iQ) - look quite similar, yet the change 
in viscosity is about four orders of magnitude p3|] . The 
other two profiles, however, greatly differ in shape, yet 
the increase in viscosity as compared to PS(18k) is only 
two (four) orders of magnitude for PS(lOlk) (PS(600k)). 
These results suggest that the alteration in profile shape 
is caused by an increase in viscoelastic properties, i.e., the 
ratio of viscosity and elastic modulus at the typical shear 
rates to which the material is subject upon dewetting. 

To corroborate this hypothesis, we performed dewet- 
ting experiments with an entangled polymer melt of a 
single chain length and varied the temperature. This is 
expected to change the viscoelastic response in a way 
completely different from changing the molecular weight. 
We prepared 50(2) nm thick PS(65k) films on type A 
wafers. The samples were heated up to different tem- 
peratures for a certain time, namely until the holes had 
reached a radius of about 5 /im. That way, we changed 
the viscosity from rj^ A-IO'^ Pa-s (T=145°C) to 77= 5-10^ 
Pa-s (T=125°C). The profiles are shown in the inset of 
Fig. 4. We have normalized the data to account for the 
slightly different hole sizes. As can be seen in the inset, 
the profiles are similar, as opposed to the large figure. 

Since the substrate is identical for all polymer films 
shown, the friction at the interface between melt and 
substrate remains unchanged for all polymer films. The 
only difference is that with increasing chain length the 
viscoelasticity of the polymer melt plays a more and more 
pronounced role [|4|. In a theoretical calculation of the 
stability condition of thin polymeric films, Safran and 
Klein |^ found that viscoelasticity stabilizes a liquid 
film against surface undulations, which corroborates our 
results. 

In summary we may say that we report for the first 
time on an undulatory behavior of the rim of a dewetting 
polymer film. Most experiments done before with dewet- 
ting films were performed on polymers of much higher 
molecular weight than 2 kg/mol. Hence, the viscoelastic 
properties of the melt have suppressed the undulation. 
We introduced a simple model that is able to explain the 
qualitative behavior of rim and trough of a short-chain 
polymers melt. A quantitative comparison of theoretical 
and experimental results reveals a difference of roughly a 
factor of 2. For long-chain polymers, the model fails com- 
pletely due to the absence of the trough. Analyzing the 



reason for the factor 2 or for the complete failure of the 
model, respectively, we came to the conclusion that the 
viscoelastic properties of the polymer melt are to blame. 
In is interesting to note here that on the rim profiles 
an influence of the elastic properties of the melt can be 
found for chain lengths at which bulk rheometers detect 



purely viscous properties |24|. The rim profile therefore 
can be viewed as an extremely sensitive rheometer for 
small volumes on nm scales. 
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FIG. 1. Inset: AFM scan of a hole in a 6.6(2) nm thick 
PS(2.24k) film on a wafer B, 30 min at T = 80 °C. The large 
diagram shows a radial cross section of this hole, where A 
describes the height of the rim, B the depth of depression 
and L the length from the maximum of the rim to the point 
where the rim crosses the level of the unperturbed film. 

FIG. 2. a) Cross sections of a growing hole in a 11.0(4) nm 
thick PS(2.05k) film on a type B wafer (191.2(2) nm thick 
oxide layer) at a temperature of T = 90° C, as gained from an 
in situ AFM scan. Note that the profiles have been shifted 
laterally such that the three phase contact 'points' fall onto 
one single point. The maximum slope at the level of the 
unperturbed film thickness is tana = 0.026. b) Depth of 
depression B over height of rim A of the hole and c) 'wet' 
contact angle a plotted over height of rim A of the hole. 

FIG. 3. Ratio (B/A) over rim height A. Dots: experi- 
mental data for a 11.0(4) nm thick PS(2.05k) film on a type 
B wafer at a temperature of T=90°C. Lines: theoretical ex- 
pectation in lubrication approximation for the same system 
gained from the analytic calculation without any elastic prop- 
erties. Dotted line: solution for the early stage B/A 0.163. 
Dashed line: solution for the later stage using reduced vis- 
cosity rj — 300 Pas with the boundary condition gained from 
experiment: atheo = 1.0°. 

FIG. 4. Cross sections as gained from AFM scans of holes 
in polystyrene films with different chain length on wafers of 
type B: solid line: PS(2k), dashed line: PS(18k), dot-dashed 
line: PS(lOlk), dotted line: PS(600k). The thickness /i of the 
films and the diameters of the holes are chosen to be about the 
same. Inset: Hole profiles of PS (65k) taken a temperatures 
between 125 - 145 °C. Differences in film thickness and hole 
size lead to slight differences of the profiles. 
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